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Abstract. This is the fourth paper in a series studying star formation rates, stellar components, metallicities, 
and star formation histories of a blue compact galaxy (BCG) sample. Using Ha, [On]A3727, infrared (IR), radio 
(1.4 GHz) luminosities and neutral hydrogen (Hi) gas masses, we estimated star formation rates and gas depletion 
timescales of 72 star-forming BCGs. The star formation rates of the BCGs in our sample span nearly four orders 
of magnitude, from approximately 10 -2 to 10 2 Moyr -1 , with a median star formation rate of about 3M0yr _1 . 
The typical gas depletion timescale of BCGs is about one billion years. Star formation could be sustained at the 
current level only on a timescale significantly lower than the age of the universe before their neutral gas reservoir 
is completely depleted. 

To assess the possible systematic differences among different star formation rate indicators, we compared the 
star formation rates derived from Ha, [Oll]A3727, IR, and radio luminosities, and investigated the effects from 
underlying stellar absorption and dust extinction. We found that subtracting underlying stellar absorption is 
very important to calculate both dust extinction and star formation rate of galaxies. Otherwise, the intrinsic 
extinction will be overestimated, the star formation rates derived from [O n] A3727 and Ha will be underestimated 
(if the underlying stellar absorption and the internal extinction were not corrected from the observed luminosity) 
or overestimated (if an overestimated internal extinction were used for extinction correction). After both the 
underlying stellar absorption and the dust extinction were corrected, a remarkably good correlation emerges 
among Ha, [On]A3727, IR and radio star formation rate indicators. Finally, we find a good correlation between 
the measured star formation rate and the absolute blue magnitude, metallicity, interstellar extinction of BCGs. 
Our results indicate that faint, low-mass BCGs have lower star formation rates. 
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1. Introduction 

Understanding the star formation history of the universe 
is the primary goal of much current research in astron- 
omy (Madau etal. 1996; Steidel etal. 1999). Star forma- 
tion rate (SFR) is a crucial ingredient to understand the 
star formation history of galaxies at all redshifts, as well 
as the global evolution of the Universe. To obtain this un- 
derstanding, a reliable estimate of the SFR in individual 
galaxies is required. Many calibrations of SFR depend on 
the luminosity measured at various wavelengths, includ- 
ing radio, infrared (IR), ultraviolet (UV), optical spectral 
lines (such as Ha, [On]A3727) and continuum (Kennicutt 
1998). Using these SFR indicators, a number of studies 
of the SFR were performed at different redshift regimes 
(such as Gallego etal. 1995; Pettini etal. 1998; Sullivan 
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etal. 2001). Unfortunately, the agreement among these 
SFR indicators at different wavelengths is poor; the under- 
lying reasons for the differences are not well understood 
(Chariot etal. 2002). 

Blue compact galaxies (BCGs) are characterized by 
their blue color, compact appearance, high gas content, 
strong nebular emission lines , and low chemical abun- 
dances jKunth fc Ostlin 2000|l . In recent years BCGs have 
attracted a great deal of interest and have become key in 
understanding fundamental astrophysical problems. Most 
of the work carried out has focused on a statistical 
analysis of BCG samples by means of surface photom- 
etry (Papaderos et al. 1996; Vitores etal. 1996; Doublicr 
et al. 1997; Ostlin et al. 1999; Cairos et al. 2001; Gil de Paz 
etal. 2003). The measured broadband colors are strongly 
affected by interstellar extinction and gaseous emission; 
therefore, optical broadband photometry alone does not 
allow us to derive their physical properties accurately. 
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Spectroscopic information is required since spectra of 
galaxies contain a wealth of additional information con- 
cerning the physical properties of galaxies. The contin- 
uum and the absorption lines provide information con- 
cerning the stellar content, while the nebular emission 
lines provide a measure of the SFR and the interstellar 
metallicity. During the last three decades, a number of 
major objective-prism surveys for extragalactic emission- 
line objects have been carried out. Examples include 
the University of Michigan (UM) survey (MacAlpine, 
Smith, & Lewis 1977), the Universidad Complutense de 
Madrid (UCM) survey (Zamorano etal. 1994, 1996), the 
KPNO International Spectroscopic Survey (KISS) (Salzer 
etal. 2000). Some BCG candidates and other emission 
line objects have been discovered. Since the objective- 
prism spectra have very low spectral dispersion and small 
wavelength range, they cannot be used to investigate 
the physical properties of BCGs. Spectroscopic follow- 
up observations were carried out to study the character- 
istics of the survey BCG candidates IjSalzer et al. 1 989 
ITerlevich et al. 19911 |Gallego etal. 1996| ). The observa- 
tions were obtained with different telescopes at different 
observatories with different equipment and instrument se- 
tups. The spectra from these observations have relatively 
low dispersions and low S/N ratios. 

We have undertaken an extensive long-slit spectral ob- 
servation of 200 BCGs since 1997. Spectra were obtained 
at the Cassegrain focus of the 2.16 m telescope of National 
Astronomical Observatory of China. A 300 line mm -1 
grating was used, giving a dispersion of 4.8 Apixel -1 . All 
200 spectral observations were finished by February 2004, 
and 97 luminous BCG spectra were studied in this series. 
The observations have been performed with the same in- 
struments; the reductions and analyses are homogeneous. 

The blue color, high gas content and strong nebular 
emission-line spectrum of BCGs indicate that they have 
a high star formation activity. One of the key parameters 
in developing an understanding of the formation and evo- 
lution of BCGs is the determination of the SFRs. While 
SFRs are central to discussions of the evolution of BCGs, 
there are surprisingly few quantitative SFRs for BCGs 
available in the literature. In addition, the values avail- 
able in the literature suggest that the SFRs of BCGs 
might vary over five orders of magnitude, and the mean 
SFRs given in the existing investigations are noted to also 
vary over two orders of magnitude. It is therefore not clear 
whether and to what extent the wide range of BCG SFRs 
results from the use of different techniques, assumptions, 
or sample selections, employed in the available studies. 

Furthermore, the studies of intermediate red- 
shift compact galaxies in the Hubble deep field 
HGuzman et al. 19961 |Phillips et al. 1997| ) reveal that 
these include galaxies with properties very similar to 
those of local luminous BCGs, such as small sizes, high 
luminosities and blue colors. Because spectroscopic 
observation for these intermediate redshift galaxies is 
difficult (for example, Ha is shifted outside the optical 
range), some other SFR indicators or methods must be 



used to determine the SFRs of these intermediate red- 
shift galaxies IjGuzman et al. 19971 IHammer et al. 20011 
IGuzman etal. 2003J) . Therefore, estimating SFRs of 
nearby BCGs with different SFR indicators may help us 
in accurately estimating SFRs of intermediate redshift 
galaxies and studying their evolutionary history. 

Based on our homogeneous BCG optical spectral sam- 
ple, we calculate SFRs and gas depletion timescales of a 
sample of star-forming blue compact galaxies, using the 
Ha and [On]A3727 emission line fluxes and the IR and 
radio luminosities. We compare the global SFRs derived 
from [Oll]A3727 and Ha, with or without corrections for 
dust extinction and underlying stellar absorption, to the 
SFRs derived from the IR and radio luminosities of these 
galaxies and investigate the source of the discrepancies 
among these different SFR indicators. We also investi- 
gate the relationship between the intrinsic SFR and other 
galaxy properties. We found that subtracting underlying 
stellar Balmer absorptions is very important in estimating 
dust extinction and SFR of galaxies. 

The paper is organized as follows. In Sect. 2 we review 
the spectroscopic observation and reduction, and describe 
the IRAS data, the radio 1.4 GHz observations. In Sect. 
3 we calculate multiwavelength SFR estimates for the 
galaxy sample using the available flux data. In Sect. 4, we 
compare the results obtained using the different SFR esti- 
mators and investigate the underlying reason for their dis- 
crepancies. We present some discussion in Sect. 5 and our 
conclusions in Sect. 6. Throughout this paper we assume 
Ho = 71 kms -1 Mpc -1 , and adopt the Salpeter initial 
mass function with mass ranges 0.1 Mq < M < 100 Mq. 



2. The data 

2.1. Ha and [On]A3727 emission line fluxes 

To investigate the current SFRs, stellar components, 
metallicities and star formation histories of BCGs, we have 
observed an atlas of optical spectra of the central regions 
of 97 blue compact galaxies in the first paper of this series 
( |Kong fc~C hcng 2002a). Since we want to combine the op- 
tical spectra we obtained with H I data to constrain simul- 
taneously the stellar and gas contents of BCGs, our sample 
galaxies were selected from large Hi surveys (Gordon & 
Gottesman 1981; Thuan & Martin 1981). The spectra were 
obtained at the 2.16 m telescope at the XingLong Station 
of the National Astronomical Observatory of China. A 300 
line mm -1 grating was used to achieve coverage in the 
wavelength region from 3580 to 7400 A with about 4.8 A 
per pixel resolution. The emission line equivalent widths 
and fluxes for our BCG sample were provided in the sec- 
ond paper of this series ( |Kong et al. 20 02b). The typical 
uncertainties of the measurements are less than 10% for 
Ha6563, H/34861, and [On]A3727 emission lines. 
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Fig. 1. Left: Rest-frame spectrum of Haro 35 on which a synthetic spectrum (thick line) from the empirical population 
synthesis method. Right: Effect of underlying stellar absorption to dust extinction ry . Ty or and Ty on were calculated 
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2.1.1. Stellar absorption and extinction correction 

The underlying stellar absorption affects accurate mea- 
surement of H-Balmcr emission line, and is crucial to con- 
strain the attenuation by dust and the SFRs of galax- 
ies. Some previous works have discussed this effect on ex- 
tinction and SFR estimation. Rosa-Gonzalez et al. (2002) 
found the SFRh q is close to the SFRjr while both SFR[ n] 
and SFRtjv show a clear excess for 31 nearby star- 
forming galaxies. The underlying B aimer absorption re- 
sults in an underestimate of the emitted fluxes and an 
overestimation of the internal extinction of the galaxy. 
Subtracting the underlying absorption from emission lines 
(with EW afcs (Ha) = EW afcs (H/3) = EW abs (H 7 )) brings 
into agreement the SFR in the optical with those in 
the FIR. Using 149 galaxies from the Stromlo-APM red- 
shift survey and using Chariot & Longhetti (2001) model 
(which provides a physically consistent description of the 
effects of stars, gas and dust on the integrated spectra 
of galaxies), Chariot etal. (2002) found the contamina- 
tion of Balmer emission by stellar H-Balmer absorption is 
important in estimating SFRs of galaxies with small ob- 
served Ha emission equivalent widths. Flores etal. (2004) 
studied SFRs for a sample of 16 distant galaxies; the un- 
derlying Balmer absorption lines were estimated by the 
Bica & Allion (1986) method. After subtracting the un- 
derlying Balmer absorption line effect, the extinction es- 
timates using the H/3/H7 and the Ha/H/3 and the SFR 
estimates (for galaxies with SFRir below ~100 M Q yr _1 ) 
using SFRhq and SFRir are in excellent agreement. 

To derive the equivalent widths for these underlying 
stellar absorption lines and correct the measured Balmer 
emission line fluxes for these absorptions, we have ap- 
plied an empirical population synthesis method to our 



BCG spectra l|Cid Fernandes etal. 2001|> . The left panel 
of Figure ^ illustrates the spectral fit of one of our sam- 
ple galaxies. The synthetic spectrum (thick line) from 
an empirical population synthesis analysis provides good 
fits to the observed spectrum of BCGs (thin line). The 
absorption wings of H/3 and H7 in the observed spec- 
trum is well reproduced by the synthetic spectrum too 
( |Kong et al. 2003"| ). The synthetic stellar spectrum is then 
subtracted from the observed spectrum to yield a pure 
emission spectrum uncontaminated by underlying stel- 
lar absorption, which we can use to measure the Balmer 
decrement, derive extinction and estimate SFR. 

In the right panel of Figure ^ we show the effect of 
underlying stellar absorption on the internal extinction 
of galaxies (ry = 0.921Ay is the effective absorption 
optical depth at 5500A by interstellar dust), where Ty 
was determined using these two strongest Balmer lines, 
Ha, H/3, and the Chariot & Fall (2000) effective absorp- 
tion curve t\ — Ty (A/5500A)~ 07 (which is similar to 
the Calzetti (2001) effective absorption curve at optical 
wavelengths; fc(Ha) - fc(H/3) = -0.57 for the Chariot & 
Fall (2000) effective absorption curve is almost same as 
fc(Ha) - fc(H/3) = -0.58 for the Calzetti (2001) effective 
absorption curve and varying the extinction curves has 
negligible effects in the visible range) . T v on was calculated 
from the Fka and Fk/3, where the underlying stellar ab- 
sorption was not subtracted. T v or was calculated from the 
Fh q and Fh/3, and the underlying stellar population ab- 
sorption was subtracted. From this figure, we find that 
subtracting underlying stellar absorption is very impor- 
tant in calculating the Ty; otherwise the internal extinc- 
tion of galaxy will be overestimated. Using t, 
nal extinction of each galaxy was corrected. 
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2.1.2. Aperture corrections 

In addition to the underlying absorption and extinction 
correction, the emission line fluxes also require an aper- 
ture correction to account for the fact that only a limited 
amount of emission from a galaxy is detected using a slit 
2"~ 3" wide and 3' long. 

Two different methods were used to estimate the 
slit covering fraction. The first method is similar to 
that used by Kewley etal. (2002). We convolve the 
galaxy spectral energy distribution with the normalized 
B-filter response curve (Bessell 1979) in £?-band magni- 
tude, mB(spec), and estimate the fraction, /, by compar- 
ison with the B-band photometry for the whole galaxy: 
f = l0 -0A[ mB -m B ( S pec)] mB ig the B _ band ma g n itude 

for the whole galaxy, which is taken from Gordon & 
Gottesman (1981), de Vaucouleurs etal. (1991, RC3) and 
the NASA/IPAC Extragalactic Database (NED). 

The second method is similar to that used by Carter 
etal. (2001). We use the digitized Palomar (POSS) im- 
ages of BCGs to estimate the slit covering fraction for 
the galaxies. We measure the ratio of the photographic 
flux within the slit aperture to the total photographic 
flux, / = Fsut/Ftotai, where F s i it is the photographic 
flux within the slit aperture, and F to tai is the total photo- 
graphic flux of a galaxy. 

The results from these two different methods are in 
good agreement in most cases. We use the average slit 
covering fraction from these two methods for the aperture 
correction. Most of our spectra contain 20—50% of the to- 
tal flux of the galaxies; for the largest galaxies this fraction 
can be as small as 10%; for the smallest galaxies this frac- 
tion can be as large as 90%. The average slit fraction of our 
sample is 37%. The absolute fluxes of Ha and [O n]A3727, 
F a bs, were calibrated from the spectral flux, F spec , and the 
average slit covering fraction with F a b s = F spec / f. 

2.2. Infrared data 

Indicators of SFR at longer wavelengths are less effected 
by dust extinction than those at UV or optical wave- 
lengths. To investigate the underlying reasons for the dis- 
crepancies between different SFR indicators, we also use 
the infrared and radio luminosities to derive the SFR of 
BCGs. 

Infrared Astronomical Satellite (IRAS) products in- 
clude 12, 25, 60 and 100 /im flux densities for galax- 
ies from the Point Source Catalog (brighter than about 
0.5 Jy) and Faint Source Catalog (typically brighter than 
0.2 Jy) (Moshir et al. 1992). We cross-correlated our galax- 
ies with the IRAS Faint Source Catalogue and Point 
Source Catalogue. The IRAS beam size at 60 [im is 1.5', 
and the IRAS positional uncertainty is 30". Because the 
spectroscopic positional uncertainty is only a few arcsec- 
onds and because all BCGs are smaller than 2' in extent, 
we used a detection radius of 30". Any IRAS source asso- 
ciated with a BCG should be detected within this radius. 
We found that of the 72 BCGs in our sample, 58 are de- 



tected at all four bands by IRAS within 30"; most of them 
have moderate or good quality fluxes at 60 and 100 /um. 

For each galaxy that was detected by IRAS, the total 
far-infrared flux Fir = F(8 — 1000 fim) was estimated from 
the observed IRAS flux densities at 25, 60, and 100 /im, 
using the method of Dale & Helou (2002). The total IR 
flux determined in this way is typically 2.2 ± 0.4 times 
larger than the FIR estimator of Helou et al. (1985), and 
similar to that from the estimator of Sanders & Mirabel 
(1996). 

2.3. 1.4 GHz radio flux 

The insensitivity of radio wavelengths to dust obscura- 
tion makes radio emission a particularly attractive way 
of estimating SFRs in star-forming galaxies. To calculate 
the SFRs of BCGs, we extract the radio 1.4 GHz inte- 
grated flux for our sample galaxies from the NRAO VLA 
Sky Survey (NVSS) and Faint Images of the Radio Sky at 
Twenty-cm (FIRST) survey. 

The NVSS is a 1.4 GHz continuum survey covering the 
entire sky north of -40 deg declination, made using the 
Very Large Array (VLA) operated by the National Radio 
Astronomy Observatory. About 220,000 individual snap- 
shots (phase centers) have been observed. The resulting 
images have a rms noise level of about 0.5 mJy beam -1 
and a FWHM resolution of 45". A detailed description is 
given in Condon et al. (1998). For 45 galaxies in our BCG 
sample, an NVSS catalog source is close enough (< 30") 
to be a probable identification. 

The FIRST is a 1.4 GHz radio continuum survey cur- 
rently covering about 9033 square degrees of the North 
Galactic Cap. The images have a typical rms noise level 
of 0.15 mJy with a FWHM resolution of 5" (White 
etal. 1997). From the FIRST catalog (03aprll version, 
2003 April 11), 24 galaxies were matched with our BCG 
sample (< 30"). For these 24 FIRST galaxies, 22 of them 
were observed by the NVSS too. Because of the much 
higher resolution of FIRST compared to NVSS, we will use 
the FIRST flux density for SFR calculation when galaxies 
were observed by both FIRST and NVSS. 

2.4. Distances 

Of the 72 star-forming BCGs in our sample, all have 
[O n]A3727 and Ha fluxes, 58 have IR fluxes, and 47 have 
1.4 GHz radio fluxes. To calculate the SFRs of galaxies, we 
convert the fluxes into luminosities, using L = 4ttD 2 x F. 
Here D is the cosmological luminosity distance, calcu- 
lated using the heliocentric velocities of the galaxies and 
a Hubble constant value of Ho = 71kms~ 1 Mpc~ 1 from 
WMAP (Spergel etal. 2003). 

3. Multiwavelength star formation rates 

In this section, we calculate the SFRs of 72 BCGs from 
the luminosities of [On]A3727 and Ha nebular lines, IR 
and radio continuum using well-known calibrations. 
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Fig. 2. Histograms showing the range of derived SFRs 
for the star- forming BCGs. It indicates the broad range 
of SFRs exhibited in blue compact galaxies, from approx- 
imately 10 -2 to several times lOMQyr -1 . The median 
SFRs from [On]A3727, Ha, IR and radio luminosities are 
2.8, 2.4, 3.6 and 5.3M Q yr _1 , respectively. 

3.1. SFR from Ha luminosities 

Hydrogen recombination line fluxes have been used very 
extensively to estimate the SFR, since they are propor- 
tional to the number of photons produced by the hot stars, 
which is in turn proportional to their birthrate. Most ap- 
plications of this method have been based on measure- 
ments of the Ha line. Other hydrogen recombination lines, 
such as H/3 and H7, were not widely used because the ef- 
fect of the underlying stellar absorption and dust extinc- 
tion are more important in these blue lines. 

For the estimate of the SFR from Ha luminosities we 
used the expression given by Kennicutt (1998): 

SFR Ha (M q yr" 1 ) = 7.9 x 1(T 42 Lh q (ergs^ 1 ). (1) 

which is valid for a T e = 10 4 K and Case B recombination, 
i.e. all the ionizing photons are processed by the gas. 

Figure 2a shows that the SFRs derived from Ha fluxes 
for the 72 star-forming blue compact galaxies in our sam- 
ple. SFRh q of BCGs cover almost four orders of magni- 
tude, from 0.01 to 65 M Q yr -1 and with a median SFR of 
2.4M Q yr- 1 . 

3.2. SFR from [0\\]\3727 luminosities 

If recombination lines are unavailable, the collisionally ex- 
cited lines of heavy elements may be used in their stead. 
The [On]A3727 doublet is the most natural star forma- 
tion tracer for objects with redshifts larger than 0.4 where 
Ha is shifted outside the optical range. The [On]A3727 
luminosity is not directly coupled to the ionization rate, 
and their excitation is known to be sensitive to the abun- 
dance and the ionization state of the gas. However, it is 
claimed that the excitation of [Oil]A3727 is sufficiently 
stable in observed galaxies that it can be calibrated em- 



pirically with Ha as a quantitative SFR tracer (Gallagher 
etal. 1989; Kennicutt 1992). 

Figure 8a in Kong et al. (2002) shows a plot of the 
[On]A3727 emission line flux versus the Ha emission line 
flux corrected for stellar absorption in our star-forming 
BCG sample. There is a good correlation between Fjo 11] 
and F Hq . From the data in Fig. 8a of Kong etal. (2002), 
we derive a least-squares relationship of F[on]/Fna == 0.85. 
The F[oii]/Fn a of BCGs (local galaxies) in good agree- 
ment with the values of high redshift galaxies, where 
F[on]/Fn a ~ 0.9 (Hippelein etal. 2003). Adopting the 
calibration between SFR and Ha luminosity given in 
Kennicutt (1998), the SFR and Lroni can be expressed 
as: 

SFR [OII] (M yi- 1 ) = 9.3 x KT^Lpii] (ergs" 1 ). (2) 

This formula is similar to that given by Gallagher 
etal. (1989) using a sample of 75 blue galaxies, where 
SFR [0 ii] (Moyr- 1 ) = 1.0 x 10- 41 L [O n] (ergs" 1 ) for blue 
galaxies. Using the Eq. @, we calculate the SFRr OII ] for 
72 star-forming BCGs in our sample, and show the re- 
sults in Figure^. The range of SFR [OII] from 0.01 to 67 
Mq yr -1 and with a median SFR of 2.8 M Q yr -1 . 

3.3. SFR from far-infrared continuum 

A significant fraction of the bolometric luminosity of a 
galaxy is absorbed by interstellar dust and re-emitted in 
the thermal IR, at wavelengths of roughly 10-300 /im. The 
absorption cross section of the dust is strongly peaked in 
the ultraviolet, so in principle the FIR emission can be a 
sensitive tracer of the young stellar population and SFR. 
Using theoretical stellar flux distributions and evolution- 
ary models, we can then derive the SFRs of galaxies. 

The SFR of a galaxy in our sample was estimated from 
its integrated IR luminosity (Kennicutt 1998) by: 

SFR m (M yr- 1 ) =4.5 x lO^pm (erg s^ 1 ), (3) 

where Lfir refers to the luminosity integrated over the full 
mid-infrared to submillimeter spectrum (8-1000 fim). The 
SFRir for 58 BCGs in our sample with IRAS detection 
were plotted in Figure range from 0.02 to 67 Mq yr -1 
and with a median SFR of 3.6 M yr" 1 . 

3.4. SFR from 1.4 GHz continuum 

In theory, the radio emission results predominantly from 
the nonthermal supernova remnant flux, and therefore 
should be very model dependent and critically influenced 
by the density of, and magnetic field in, the interstellar 
gas. However, the FIR/radio (the 1.4 GHz radio contin- 
uum flux) correlation suggests that the constant of pro- 
portionality is nearly the same for all star forming galaxies 
and is insensitive to unknown variables such as magnetic 
field strength (Condon 1992; Yun etal. 2001). The real- 
ity of this correlation has been repeatedly demonstrated 
by observation, not only in terms of the integrated flux, 
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but also in terms of spatial correlations between these two 
quantities within individual galaxies. Therefore, the radio 
continuum flux can be used as a star formation indicator. 

For an individual star-forming galaxy, the SFR is di- 
rectly proportional to its radio luminosity at 1.4 GHz 
(Condon 1992): 

SFR1.4 (M yr 1 ) = nx2.2x KT 22 Li. 4 (W Hz" 1 ). (4) 

Condon (1992) derives this relation by calculating the 
synchrotron radio emission from supernova remnants and 
the thermal radio emission from HII regions. Both the 
thermal and non-thermal components of the radio expres- 
sion are proportional to the formation rate of high-mass 
stars (M > 5M©) which produce supernova and large HII 
regions, so the factor n is included to account for the mass 
of all stars in the interval 0.1 — 100M©. We have assumed 
throughout a Salpeter IMF (x = 2.35), for which n = 5.5. 

The SFR1.4 for 47 BCGs in our sample with NVSS 
and FIRST detection were plotted in Figure F2Ji, range 
from 0.01 to 129 M yr" 1 and with a median SFR of 5.3 
Mq yr _1 . Most of them are consistent with those derived 
from Ha, [On]A3727 and IR luminosities, but the SFR1.4 
of I Zw 56 is much larger than those derived from other 
SFR indicators. 

3.5. SFR of BCGs 

Based on the FIR (60 fj,m) and radio continuum (1.4 GHz) 
luminosities, Izotova et al. (2000) derived the SFRs of 27 
BCGs, with a mean of about lSMoyr" 1 . Using spectro- 
scopically derived H^ luminosities, Popescu etal. (1999) 
determined the SFRs of a sample of BCGs, with a mean 
SFR of about 0.5 M© yr" 1 . Therefore, the values available 
in the literature suggest that the mean SFRs given in the 
existing investigations vary over two orders of magnitude. 
It is therefore not clear whether and to what extent the 
wide range of BCG SFRs results from the use of different 
techniques, or sample selections, employed in the available 
studies. 

In this section, we have determined the SFRs of 72 blue 
compact galaxies, using the Ha, [On]A3727 emission line 
fluxes, the IRAS, and the 1.4 GHz radio continuum lumi- 
nosities. Since we have determined the SFRs for our ho- 
mogeneous BCG sample, and calculated the median SFR 
of the sample using these different SFR indicators, it can 
be used to understand the reasons for these discrepancies 
in the literature. 

Fig. [21 shows that the SFRs derived from different star 
formation indicators cover almost four orders of magni- 
tude, from 10~ 2 to 10 2 M©yr _1 . After correcting for the 
underlying stellar absorption and the dust extinction ef- 
fects, the SFR of each galaxy derived from those different 
SFR indicators is similar. Using different techniques will 
not result a large different median SFR of galaxies if we 
have observed the same sample galaxies over the entire 
wavelength bands. 



From Fig. FJ1 we also found that the medium and av- 
erage SFRs from these different indicators of our sam- 
ple galaxies are different; values derived from the IR and 
radio data are much larger than those derived from the 
optical emission lines. Therefore, the sample selection ef- 
fect is an important reason for the different mean SFR of 
BCGs in the literature. This can be explained by: the posi- 
tional accuracy of far-infrared observation being not good, 
therefore, it is difficult to identify an object with weak IR 
emission. On the opposite side, radio observations typi- 
cally have better positional accuracy, but the detector in 
the radio band is not sensitive. Therefore both far-infrared 
and radio observations are difficult to use to identify an 
objects with weak star formation activity; these survey 
are biased towards the galaxies with high SFRs, and the 
mean values of these survey samples are larger than those 
of optical observations. 

4. Comparison of SFRs derived using different 
standard estimators 

We have applied the commonly used SFR estimators to 
the 72 star-forming BCGs in our sample. In this section, 
we will investigate the SFR relation from these different 
standard estimators and the underlying reason for their 
discrepancies. 

4.1. SFR [ou] and SFR Ho 

Figure |3| shows the logarithmic difference between SFRjj q 
and SFR[on] for the 72 BCG in our sample as a function 
of the B-band absolute magnitude, Mb- In each panel, we 
indicate the median A log(SFR) offset and the rms scatter. 

In Fig. Et, SFR[ ii] and SFRjja were calculated with 
equations (1), (2) and the observed £[011] and Lh q data, 
where the underlying stellar absorption and the internal 
extinction were not corrected. Fig. Et is consistent with 
the recent result by Chariot et al. (2002) that bright galax- 
ies tend to have lower [On] /Ha ratios than less luminous 
galaxies. We find that, for M B < -18 mag, [On]A3727- 
derived SFR estimates are typically 30 per cent smaller 
than Ha-derived ones, while at fainter magnitudes, the 
two estimators give similar results. At fixed magnitude, 
the rms scatter in SFR[on]/SFRH a is a factor of ~ 1.7. 
The scatter in this plot and the weak trend with magni- 
tude can probably be explained by variations in the effec- 
tive gas parameters (ionization, metallicity, dust content) 
of the galaxies. 

In Fig. Et>, SFRs were calculated from the L[oii] an d 
Luct data for which the internal extinction was corrected 
by Ty on , which was determined by the observed Ha/H/3 
of the galaxy (the underlying stellar absorption was not 
subtracted). In star-forming galaxies, the Balmer emis- 
sion lines from the ionized gas appear superimposed on 
the stellar absorption lines; this effect grows in impor- 
tance towards the higher order Balmer lines, and the re- 
sult is an underestimate of the emitted fluxes and an over- 
estimate of the internal extinction, as show in Fig. 1. 
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Fig. 3. Logarithmic difference between SFRna and 
SFR[ OII ] plotted against absolute rest-frame Mb magni- 
tude for star-forming BCGs. SFRs derived from Ln a and 
L[ou] : ( a ) before applying any correction for Lh q and 
L[on]i (b) after correction for dust extinction, using the 
Ty on (see Fig.l); (c) after further correction for the un- 
derlying stellar absorption. In each panel, the median 
Alog(SFR) offset is indicated along with the associated 
mean offset and the rms scatter. 

This overestimated internal extinction correction makes 
the [O n] A3727-derived SFR estimates typically a factor 
of ^ 2 higher than Ha-derived ones. 

In Fig. the SFR[ n] and SFRh q were determined 
by the Lroii] and Lrqi where the underlying stellar ab- 
sorption was subtracted and the internal extinction was 
corrected. The SFR of BCGs from Fjon] and Fh q have 
an excellent match, which can be understood easily since 
we derived the SFRpn] formula (equation 2) from the 
relation between Fjon] and -Fh q . In Fig. Ok the rms scat- 
ter in SFR[oii]/SFRhq becomes smaller, about 1.4. The 3 
low luminosity BCGs are unaffected by stellar absorption 
and extinction. This is because: 1) They are low metallic- 
ity galaxies, and have a small dust component. The effect 
of internal extinction is not important for them. 2) They 
appear dominated by a recent burst of star formation, 
which causes their extremely strong emission line spec- 
trum. The underlying stellar absorption is much weaker 
than the emission component. In addition, the correlation 
between Fu_ a and Fjo n] derived from the whole sample 
(most of them are high luminosity galaxies) is not good 
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Fig. 4. Logarithmic difference between SFRh q and 
SFRir, plotted against absolute rest-frame Mb magnitude 
for 58 BCGs, which were observed by IRAS. SFRs derived 
from Lhq and Lfir: (a) before applying any correction 
for Lhq (b) after correction for dust extinction, using the 
Ty on , (c) after further correction for the underlying stellar 
absorption. In each panel, the median Alog(SFR) offset 
is indicated along with the associated mean offset and the 
rms scatter. 

for these low metallicity galaxies, and may underestimate 
the SFR[on] for low luminosity galaxies. 

4.2. SFR IR and SFR Ha 

Figure [1] shows the logarithmic difference between SFRh q 
and SFRir for 58 BCGs (IRAS data are available) in our 
sample as a function of Mb . 

In Fig. HJi, SFRir was derived from equation (3), 
SFRhq was calculated with equation (1) and the observed 
Ha luminosity data, Lho- The underlying stellar absorp- 
tion and the internal extinction were not corrected in the 
observed Ha luminosity data in this panel. We find the dif- 
ference between the SFRs estimated using Lhq and those 
estimated using the total FIR luminosity is large. The 
SFRs derived using Ha are lower than those obtained us- 
ing Fir. The median discrepancy tends to increase with 
B-band absolute luminosity and amounts to a factor of 2 
at Mb < — 18. The discrepancy and the scatter is similar 
to the SFRs derived using [On]A3727. The median and 
the standard deviation of the logarithmic difference be- 
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tween SFRh q and SFRir for 58 BCGs in our sample is 
—0.32 ± 0.25, which similar to (or a little smaller than) 
that of Chariot etal. (2002) (-0.47 ± 0.38) and of Rosa- 
Gonzalez etal. (2002) (-0.53 ± 0.52). 

In Fig. 2t>, SFR Ha was calculated from the L^ a data 
with the internal extinction corrected by T v on (see Fig. 1 
for more detail). Since the effect of the underlying stellar 
absorption was not considered, the median of the Ty on is 
about 1.9, which is much larger than the median of the 
intrinsic internal extinction value (r™ r ~ 0.8) of galax- 
ies. The internal extinction value of galaxies was overes- 
timated; it makes the differences between the SFRs esti- 
mated using optical lines and those estimated using the to- 
tal IR luminosities greater. Fig.0J) shows that Ha-derived 
SFR estimates are typically a factor of ~ 3 higher than 
FiR-derived ones. 

For the galaxies in our sample, the mean of the 
equivalent width of H-Balmer absorption corresponds to 
EW£ b a s w -2 A EWf$ « -4 A (Kong etal. 2003), 
which is consistent with the results from Bruzual & 
Chariot (2003), EW$£ « 0.6EW^ S . The correction for 
stellar absorption is therefore quite important for galaxies 
with narrow observed Ha emission equivalent widths. This 
is shown in Fig. Et, where we corrected the Ha-derived 
SFRs of individual galaxies in our sample for both the un- 
derlying stellar absorption and the interstellar dust extinc- 
tion (using Ty r ). These corrections bring the Ha-derived 
SFRs into agreement with the Fpu^-derived SFRs. 

The median and the standard deviation of the logarith- 
mic difference between SFRh q and SFRir for 58 BCGs in 
our sample is 0.06 ± 0.13. Our results agree with those of 
Rosa-Gonzalez et al. (2002), and also Chariot et al. (2002). 
Rosa-Gonzalez et al. (2002) developed a method to esti- 
mate simultaneously the intrinsic visual extinction Ay 
and the underlying Balmer absorption. The values of A v 
were applied to the emission line fluxes. They also found 
that the extinction correction which includes the effects of 
an underlying stellar Balmer absorption brings into agree- 
ment SFRhq and SFRir . 

4.3. SFR1.4 and SFR Hq 

Figure |S] shows the comparison between SFRs derived 
from 1.4 GHz (FIRST and NVSS) and the Ha luminosi- 
ties. If the underlying absorption and dust extinction 
were not corrected for, the trend between the SFR1.4 and 
SFRh q is similar to that of SFRir and SFRh q (see Fig.0^,, 
b), since the radio luminosities and the IR luminosities 
have a tight correlation. We show the estimated SFRs de- 
rived from Ha luminosities, with the underlying stellar 
absorption and the internal dust extinction corrections, 
as a function of SFR1.4 in Figure only. 

The Spearman rank correlation coefficient is r s = 0.9 
for the relation between SFR1.4 and SFRh q ; this implies 
that the correlation between SFR1.4 and SFRuct is signifi- 
cant at the 6(7 level. This correlation has been reported by 
other authors for both the case of nearby galaxies (such 



. NVSS 
. FIRST 



Fig. 5. Comparison of SFRh q with SFR1.4 from both 
NVSS (circles) and FIRST (triangles). The dashed line 
is the relation SFRh q = SFR1.4, and the dotted line is a 
linear least-squares fit to the data points. The Spearman 
Rank correlation coefficient between SFRhq with SFR1.4 
is 0.9. 



as Afonso etal. 2003; Hopkins etal. 2003) and for more 
distant galaxies (Flores etal. 2004). Another interesting 
point about this relation is that bright galaxies (high 
SFRs) tend to have higher SFRi.4/SFRh q ratios than less 
luminous galaxies. This can be seen clearly from the one- 
to-one line (dashed line) and the ordinary least squares 
fit line (dotted line). For those galaxies with SFRh q < 1 
M0yr -1 , SFRs derived from 1.4 GHz radio luminosities 
are less than those derived from Ha luminosities; the prob- 
able reason for this is that the detector in radio band is 
not sensitive, and the uncertainty for those faint galaxy 
observations is large. The deviation at high radio lumi- 
nosities, as well as the increasing scatter in the correlation, 
could result from a relatively large amount of extinction 
in those objects undergoing the most vigorous star forma- 
tion, from the absorption of Lyman photons by dust, or 
from an IMF that weights differently the high-mass stars 
mainly responsible for Ha and the lower mass stars that 
dominate the supernova numbers (Cram etal. 1998). 

In this section, after comparing with SFR[ n] and 
SFRh q ; SFRir, and SFRh^, we found that subtracting the 
underlying stellar absorption is very important to deter- 
mine the internal dust extinction, and also in estimating 
the SFRs of galaxies, especially high luminosity galaxies. 
Without the underlying stellar absorption and the inter- 
nal dust extinction correction, the SFRs estimated with IR 
and radio luminosities are larger than those derived from 
optical spectral lines (SFRir > SFR Ha > SFR[ n])- This 
result is consistent with what has been found and shown 
in Madau-type plots in recent years, where the SFRs ob- 
tained from UV and optical data are much lower than that 
obtained from mm and submm observations at interme- 
diate and high redshifts. In order to reach agreement be- 
tween both determinations, fixed (and somehow arbitrary) 
amounts of extinction have been applied to the UV/optical 
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late the gas depletion time-scale of BCGs, obtained by 
dividing the available gas mass by the present SFR: 
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Fig. 6. The gas depletion time-scale of BCGs from Ha, 
[O n] A3727 emission lines, far-infrared, and radio luminosi- 
ties. The typical gas depletion time-scale of BCGs is about 
a few billion years. 



data, because at the moment, the intermediate and high 
redshift samples do not allow a reliable determination of 
the dust extinction (Rosa-Gonzalez etal. 2002). If the un- 
derlying Balmer absorption was not subtracted, the dust 
extinction from Balmer lines (such as Ha, H/3) would be 
overestimated. The SFRh q and SFR[ n] show a clear ex- 
cess and the excess is much larger for SFRpii] than for 
SFR Ha (SFR m < SFR Hq < SFR [OI t]). After taking into 
account the corrections for Balmer absorptions and dust 
extinction, the SFRs in the optical with those in the FIR 
and radio show a very good agreement. 



5. Discussion 

5.1. Gas depletion timescales 

The star formation history of BCGs is still an open ques- 
tion, BCGs arc cither genuinely young galaxies, or old 
galaxies that have resumed forming stars at a prodigious 
rate. In the third paper of this series (Kong etal. 2003), 
using a technique of empirical population synthesis based 
on a library of observed star cluster spectra, we find that 
BCGs are typically age-composite stellar systems; in ad- 
dition to young stellar populations, intermediate-age and 
old stellar populations contribute significantly to the 5870 
A continuum emission of most galaxies in our sample. 
Since all BCGs have been found to contain a population 
of old stars, our results suggest that BCGs are primarily 
old galaxies with discontinuous star formation histories. 

Since the galaxies in our sample were selected from 
some large neutral hydrogen (H i) surveys (Gordon & 
Gottesman 1981; Thuan & Martin 1981), the Hi gas mass 
(Mm) of our sample galaxies was measured by these Hi 
surveys. The SFR of our sample galaxies were measured 
in Sec. 3, using the Ha, [On]A3727 emission lines, far- 
infrared, and radio luminosities. Therefore, we can calcu- 



Tgas = Mhi/SFR, 



(5) 



where Mm is the total Hi mass content and SFR is 
the ongoing star formation rate. 

Figure shows the distribution of gas depletion 
timescales for the galaxies in our sample. Most galaxies 
are consuming their interstellar gas at a rate which can 
only be sustained for another few billion years. The me- 
dian value of the gas depletion time-scale of BCGs is less 
than or about one billion years, and the gas depletion 
time-scale of BCGs is generally much shorter than the age 
of the universe. It indicates that star formation in BCGs 
could be sustained at the current level only on a time scale 
significantly lower than l/H® (cosmological time) before 
their neutral gas reservoir will be completely depleted. 

To understand the evolutionary connection between 
blue compact star-forming galaxies observed at intermedi- 
ate redshifts (0.1 < z < 1) and nearby blue compact galax- 
ies, Pisano etal. (2001) have compared Arecibo Hi 21 cm 
spectroscopy and Keck HIRES spectroscopy of 11 nearby 
blue compact galaxies. Using the SFRh q and the Hi 
mass, they calculated the gas depletion timescale, T gas = 
Mri/SFR, of each galaxy. The SFRs of BCGs range from 
0.23—23.5 M Q yr _1 , and gas depletion timescales range 
from 0.18—7.29 Gyr. Using the optical spectroscopy, deep 
optical/near-IR photometry of 4 luminous blue compact 
galaxies, and spectral evolutionary models, Bergvall & 
Ostlin (2002) concluded that star formation in luminous 
BCGs is intense, corresponding to gas depletion timescales 
of the order of 100 Myr. Our results are consistent with 
those works: the luminous BCGs have shorter gas deple- 
tion timescale, and the typical gas depletion timescale of 
BCG is about a few billion years. 

Based on surface photometry and spectroscopic analy- 
sis, most of BCGs have been found to contain a population 
of old stars. Based on the gas depletion estimation, the 
high star formation activities of BCGs must be transient. 
Therefore, the general star formation history of BCGs is 
that these galaxies undergo a few or several short bursts of 
star formation followed by longer more quiescent periods. 

5.2. SFRs versus M B 

Figure shows the SFR against absolute B-band mag- 
nitude for BCGs; the SFRs were derived from the Ha 
(circles) and from IR (triangles) luminosities, and Mb 
were taken from Kong & Cheng (2002a). There is a 
trend between B-band magnitude and SFR of BCGs, the 
Spearman Rank correlation coefficient being r s = —0.9. 
The higher SFRs are clearly being hosted by galaxies with 
higher luminosities although there is a significant scatter 
in the trend. The least-squares fit of these data, shown by 
the dash line, yields an increased SFR towards high lumi- 
nosity galaxies, SFR = -(7.65 ±0.25) - (0.42±0.01)M B . 
This relation implies that the SFR per solar luminosity in 
BCGs is constant over a range of 7 magnitudes. 
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Fig. 7. Correlation between the star formation rates de- 
rived from the infrared, Ha luminosities and the abso- 
lute blue magnitude of galaxies. The data suggest that 
there exists a good correlation between star formation 
rates of BCGs and Mb- The dashed line shows the or- 
dinary least-squares fit for BCGs. The different symbols 
indicate SFRs derived from Ha (circles) and IR (triangles) 
for BCGs in this paper, and IR (stars) for galaxies from 
Kong etal. (2004). 

In Kong etal. (2004), the authors have compiled an 
ultraviolet-selected sample of 115 nearby, non-Seyfert 
galaxies spanning a wide range of star formation activities, 
from star burst to nearly dormant, based on observations 
with various ultraviolet satellites. To check the correlation 
between SFR and Mb for other galaxy types, we over-plot 
these galaxies in Figure as stars. IRAS flux densities 
/i/(25/xm), /j/(60/Km) and /„(100/xm) were used to esti- 
mate the total far-infrared flux Lir and SFRir. The ap- 
parent blue magnitudes from De Vaucouleurs et al. (1991) 
and the redshift from the NASA/IPAC Extragalactic 
Database (NED) were used to calculate the absolute B- 
band magnitudes. From Figured we found there is a scat- 
tered but good correlation between the Mb and the SFR 
of other galaxy types. Starburst galaxies lie substantially 
above the Mb ~ SFR relation for BCGs, and the less 
vigorously star forming galaxies are located in the region 
below this relation. Our results are consistent with the 
much large sample of ~ 10 5 SDSS galaxies studied by 
Brinchmann etal. (2004): they found a clear correlation 
between SFR and stellar mass (luminosity) over a signifi- 
cant range in log M*. 

While optical spectral line fluxes, UV, IR, and radio 
luminosity are commonly used as SFR indicators, the lu- 
minosity at B-band wavelengths is also dominated in star- 
burst galaxies by young stellar populations. In the absence 
of other observations, the B-band luminosity may thus 
be used as an SFR indicator for nearby large-area im- 
age surveys (such as SDSS) and deep image surveys (such 
as COSMOS and Subaru Deep Field , Ouchi etal. 2003). 
BCGs can be selected according to their compactness and 
blue colors from images in these surveys, and the best-fit 



relation between M b and SFR can be used to estimate the 
SFR of BCG. The SFR based on Mb is not very accuracy, 
since the scatter for the Mb and SFR relation is large. 

5.3. SFRs and dust extinction 

Star forming galaxies in the local universe have been 
shown to exhibit a correlation between obscuration and 
FIR luminosity (Wang & Heckman 1996). This implies 
that the obscuration in a galaxy is related to its SFR. This 
same general effect can also be seen in other recent stud- 
ies (such as Takagi, Arimoto & Hanami 2003). Buat et al. 
(2002) found a good relation between the extinction in 
the Ha line and the FIR luminosity for starburst galaxies. 
Hopkins etal. (2001) and Sullivan etal. (2001) report a 
positive correlation between the dust extinction traced by 
the Balmcr decrement and the SFR of the galaxies traced 
by their total FIR or Ha luminosities. 

Figure|S]shows the correlation between the dust extinc- 
tion (Ty r ) and the SFR of BCGs in our sample. Using the 
data of SFRh q and SFRir, we perform a least squares fit 
to the r v or and SFR to derive SFR = -(0.41 ± 0.09) + 
(0.73 ± 0.0.08)Tf or . The Spearman Rank correlation co- 
efficient between t v ot and SFR is 0.62. There is a clear 
trend between the dust extinction and the SFR of BCGs; 
the value of SFR increases with t v ot '. We see evidence 
for a strong correlation between SFR and metallicity (in 
the next subsection), presumably reflecting a correlation 
between metallicity and dust content. The correlation be- 
tween SFR and dust may simply reflect a more fundamen- 
tal correlation between dust and metallicity. This ques- 
tion clearly merits further investigation. The scatter in 
these trends may be related to intrinsic differences within 
galaxy populations. Because greater SFRs correspond to 
greater extinction in the dust, then the SFRs of high red- 
shift galaxies derived from either Balmer fluxes or from 
UV continuum measurements would tend to be systemat- 
ically underestimated, if the dust extinction were not cor- 
rected. Since the dispersion between the dust extinction 
and the SFR of BCGs is large, the situation is worse for 
these low extinction galaxies (most BCGs have low t v ot ). 
This relation cannot be used to estimate the amount of 
dust extinction in BCGs. 

5.4. SFRs vs. Metallicity 

Figure El shows the correlation between oxygen abun- 
dance and SFR. The oxygen abundance of BCGs was es- 
timated by the N2 = [Nn]A6584A/Ha estimator from the 
results of Denicolo, Terlevich & Terlevich (2002). These 
authors found that the N2 estimator follows a linear re- 
lation with log(0/H) that holds for a wide abundance 
range, from about l/50i/i to twice the Solar value. A 
positive trend of 12+log(0/H) with SFR Hq and SFR IR 
can be seen, in the sense that galaxies with higher global 
SFRs also show higher oxygen abundance. The Spearman 
Rank correlation coefficient for this correlation is 0.67, 
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Fig. 8. The dust extinction traced by the Balmer decre- 
ment versus the SFR. Ty r were calculated from the Fn a 
and Fn/3 fluxes where the underlying stellar population 
absorption were subtracted. 
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Fig. 9. Metallicity, 12 + log(0/H), as a function of SFR 
(SFRhq and SFRir). Galaxies with higher global SFRs 
also show higher oxygen abundance. 



an unweighted linear regression yields 12 + \og(0/H) = 

0. 2410g(5'Fi?) + 8.52. 

Zaritsky, Kennicutt & Huchra (1994) found the lumi- 
nosity (mass)-metallicity relationship exists over a large 
range of luminosity galaxy type, and metallicity. A simi- 
lar relationship was seen for BCGs in our sample (see Fig. 
9 in Kong etal. 2002). This correlation reflects the fun- 
damental role that galaxy mass plays in galactic chemi- 
cal evolution. With the luminosity and SFR correlation 
(see Figured, the most straightforward interpretation of 
the correlation between SFRs and metallicity is that more 
massive galaxies form fractionally more stars in a Hubble 
time (higher SFR) than low-mass galaxies, and then have 
higher metallicity. From Figure we also found that the 
slope of the SFR and metallicity relation for these low- 
mass galaxies increases more quickly and the slope of the 
more massive galaxies increases more slowly. It can proba- 
bly be explained as in massive galaxies, which are already 
relatively metal rich, it is more difficult to significantly 
alter 12 4- \og(0/H) than it would be in a low-mass, low- 
metallicity galaxy. 

6. Summary 

In this paper, we have determined the star formation rates 
and gas depletion timescales of 72 blue compact galaxies 
and investigated the source of the discrepancies between 
these different SFR indicators. Then we analyzed the cor- 
relation between the star formation rate and the absolute 
B-band magnitude, dust extinction and metallicity for the 
star-forming BCGs in our sample. 
Our main results are: 

1. The SFRs derived from the Ha, [On]A3727 emission 
lines, the IRAS, and the 1.4 GHz radio luminosities 
show that the SFR of BCGs cover almost four orders 
of magnitude, from 10~ 2 to 10 2 Mqvi --1 . 



2. The medium and average SFRs from these different 
indicators for our sample galaxies are different. The 
values derived from IR and radio continuum data are 
much larger than those derived from emission lines in 
the optical spectrum. The sample selection effect is 
important for the different mean SFRs of BCGs in the 
literature, determined using different SFR indicators. 

3. Subtracting the underlying stellar absorption is very 
important when calculating both the dust extinction 
and SFRfj a , SFRjonj- Otherwise, the intrinsic extinc- 
tion will be overestimated or underestimated; star for- 
mation rates derived using [On]A3727, Ha, and IR 
have large discrepancies. 

4. The gas depletion time-scale of BCGs is generally 
much shorter than the age of the universe. Star forma- 
tion in BCGs could be sustained at the current level 
only on a time scale significantly lower than the cos- 
mological time before their neutral gas reservoir will 
be completely depleted. 

5. There are some good correlations between the SFR 
and the absolute B-band magnitude, dust extinction 
and metallicity for the star-forming BCGs. The galax- 
ies with lower luminosity have lower metallicity, lower 
dust extinction and lower SFR. 
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